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LONG-TERM  GOALS 

The  work  described  here  supports  the  Office  of  Naval  Research  Departmental  Research  Initiative 
(DRI)  for  research  on  the  “Impact  of  Typhoons  on  the  Western  Pacific  Ocean”  (ITOP).  The 
theme  of  the  DRI  is  to  better  characterize  and  predict  the  ocean  boundary  layer  (OBL)  and  its 
impact  on  typhoon  (hurricane)  evolution.  This  is  one  component  of  developing  improved  pre¬ 
diction  models  for  the  coupled  atmosphere-ocean-wave  system.  Cooling  of  the  sea  surface  tem¬ 
perature  (SST)  is  a  critical  coupling  variable  influencing  atmosphere-ocean  hurricane  dynamics; 
SST  is  largely  determined  by  OBL  turbulence,  surface  wave  processes,  and  mixed  layer  entrain¬ 
ment.  Our  research  goal  is  to  model  the  strongly  forced  wind  and  wave  driven  upper  OBL  using 
turbulence  resolving  large-eddy  simulation  (LES)  with  explicit  wave  effects,  viz.,  wave-current 
interactions  and  breaking  waves  and  examine  their  impact  on  ocean  mixing  during  hurricane 
events. 


OBJECTIVES 

The  specific  research  objectives  for  ITOP  are:  (1)  conduct  process  studies  using  LES  of  the 
OBL  with  different  combinations  of  time  varying  large  scale  forcings  and  surface  wave  effects 
and  examining  their  impact  on  ocean  mixing;  (2)  evaluate  and  compare  these  LES  results  with 
predictions  obtained  using  a  1-D  column  model  of  the  OBL  based  on  the  K-Prohle  Parameteri¬ 
zation  (KPP);  and  (3)  compare  our  simulation  results  with  available  observations.  Inertial  res¬ 
onance,  storm  residence  time,  and  the  larger  scale  environment  are  some  of  the  processes  to  be 
examined  in  our  simulations. 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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APPROACH 


In  order  to  examine  the  dynamics  of  the  upper  ocean  boundary  layer  under  high  wind  forcing 
we  developed  an  LES  model  of  the  OBL  that  accounts  for  surface  wave  effects  through  non¬ 
conservative  breaking  waves  and  phase-averaged,  conservative  wave-current  interactions  which 
lead  to  Langmuir  circulations  (Sullivan  et  al.  2007).  The  equations  for  the  resolved  flow  compo¬ 
nents  are  the  Craik-Leibovich  (CL)  theory  with  crucial  wave-current  coupling  through  the  vor¬ 
tex  force.  The  larger  scale  forcing  (momentum  and  scalar  fluxes)  and  wave  fields  are  externally 
imposed  in  the  LES. 


WORK  COMPLETED 

Meetings:  We  attended  a  PI  meeting  in  Taichung  Taiwan  that  focused  on  refining  elements 
of  the  ITOP  field  campaign.  The  discussion  concentrated  on  the  deployment  of  the  numerous 
drifters,  buoys,  and  floats  and  the  coordination  with  the  aircraft. 

OBL  LES:  We  concentrated  on  adapting  our  large-eddy  simulation  code  for  the  OBL  to  a  high 
wind  regime  with  time  varying  winds,  surface  scalar  fluxes,  and  wave  fields.  For  our  LES  pro¬ 
cess  studies,  Hurricane  Frances  is  selected  as  a  canonical  storm  (e.g.,  see  Sanford  et  al.  2007). 
Hurricane  Frances  had  a  radius  of  maximum  winds  R  —  40  km  with  a  surface  rotational  wind 
speed  Vr  ~  50ms-1  (see  Zedler,  2007). 

The  general  design  of  the  LES  experiments  is  to  pass  a  hurricane  vortex  over  small  LES  do¬ 
mains  as  sketched  in  figure  1.  Then  the  surface  fluxes  and  wave  fields  imposed  at  the  water  sur¬ 
face  vary  with  time  depending  on  the  storm  propagation  speed  Vs  and  the  location  of  the  LES 
domain  ( xleSiUles )  within  the  storm  track  (see  discussion  in  RESULTS).  Our  process  studies 
examine  the  impact  of  resonant  and  anti-resonant  wind  and  inertial-current  forcing  by  locating 
LES  domains  to  the  right  and  left  of  the  storm  track. 

The  initial  temperature  state  of  the  OBL  for  all  the  LES  experiments  is  neutrally  stratified  over 
the  region  -33  m  <  z  <  0  m  bounded  by  a  stable  inversion  of  0.05  C  m-1  for  z  <  —33  m.  The 
LES  domain  (XL,  YL,  ZL)  =  (750,  750,  — 195)m  is  discretized  using  (Nx,  Ny,  Nz)  =  (500,  500, 160) 
gridpoints;  the  vertical  spacing  varies  smoothly  from  Az  =  1.0  m  near  the  water  surface  to 
about  Az  =  1.4  m  near  the  lower  boundary.  Each  simulation  is  first  run  from  a  cold  start  for 
a  period  of  ~  10  hours  to  generate  turbulence.  The  LES  is  then  restarted  using  these  archived 
turbulent  flow  fields  as  initial  conditions.  For  LES  experiments  with  wave  effects,  we  set  the 
turbulent  Langmuir  number  Lat  =  (0.3, 0.5,  oo)  and  wave  age  =  0.6  which  fixes  the  parame¬ 
ters  in  an  exponential  Stokes  drift  profile  (Sullivan  &  McWilliams,  2010).  More  than  190,000 
timesteps  are  needed  to  cover  the  entire  60  physical  hours  of  interest  for  cases  located  along  the 
path  of  maximum  wind  speed.  We  note  that  the  OBL  LES  code  used  here  is  a  variant  of  the 
new  highly  parallel  atmospheric  code  that  is  being  developed  for  the  high  resolution  air-sea  in¬ 
teraction  (HIRES)  DRI1.  The  algorithm  is  based  on  an  incompressible  Boussinesq  flow  model 
and  solves  the  governing  equations  utilizing  a  mixed  finite-difference  pseudospectral  scheme. 


1For  a  description  of  Hi-Res  see  http://airsea.ucsd.edu/hires/. 


The  parallelization  is  accomplished  utilizing  the  Message  Passing  Interface  and  a  2-D  domain 
decomposition,  see  Sullivan  &  Patton  (2008)  for  further  details. 

RESULTS 

A  typical  time  history  of  wind  and  thermal  forcing  applied  at  the  top  of  the  simulations  is  given 
in  figure  2  for  an  LES  domain  located  55  km  to  the  east  and  700  km  north  of  the  storm  cen¬ 
ter  depicted  in  figure  1.  In  this  simulation,  the  hurricane  vortex  propagates  northward  with  a 
speed  of  5.5  rns^1  and  thus  the  maximum  winds  occur  at  35.4  hours  into  the  simulation;  the 
maximum  wind  speed  is  about  43.5  ms-1.  At  this  location,  i.e.,  to  the  right  of  the  storm  track, 
the  winds  and  wind  stress  rotate  clockwise  with  a  time  scale  approaching  an  inertial  period 
T  =  27 r/f.  Inspection  of  figure  2  indicates  that  over  the  time  period  0  to  10  hours,  the  winds 
are  modest  0(5  ms-1)  while  the  surface  heat  fluxes  are  0(100  W m-2).  Also  the  wave  forcing 
is  small  because  of  small  values  of  Stokes  drift.  Hence  shear,  convection,  and  wave-current  in¬ 
teractions  all  play  a  role  in  forcing  the  OBL.  Beyond  10  hours,  shear  forcing  and  wave-current 
interactions  dominate  because  of  the  rapid  increase  in  the  surface  winds. 

The  presence  and  impact  of  Langmuir  cells  in  our  solutions  is  identified  using  flow  visualiza¬ 
tion  and  quantified  using  statistical  measures.  Vertical  velocity  w'  is  one  of  the  most  robust  and 
clear  indicators  of  Langmuir  circulations  in  simulations  with  wave-current  interaction  ( e.g see 
McWilliams,  et  al.  1997).  Measurements  of  vertical  velocity,  obtained  from  near  neutrally  buoy¬ 
ant  Lagrangian  floats  (Tseng  &  D’Asaro,  2004),  can  also  be  used  to  detect  the  presence  of  Lang¬ 
muir  circulations.  Figure  3  shows  contours  of  normalized  vertical  velocity  w' /w*  in  x  —  y  planes 
near  the  water  surface  £  ~  —5  m  at  different  simulation  times;  the  selected  times  bracket  the 
time  of  maximum  winds.  Each  panel  also  shows  a  wind  vector  oriented  in  the  mean  wind  direc¬ 
tion.  Inspection  of  the  images  shows  coherent  structures  elongated  along  the  axis  of  the  mean 
wind  direction  which  we  identify  as  Langmuir  cells.  The  overall  impression  is  that  these  coher¬ 
ent  structures  are  geometrically  similar  to  their  counterparts  found  at  much  lower  wind  speeds, 
e.g.,  see  McWilliams  et  al.  (1997).  In  the  present  solutions  the  cells  rotate  with  the  time  vary¬ 
ing  wind  direction  but  still  display  mergers  and  junctions  as  in  simulations  with  steady  forcing. 
Notice  also  that  at  the  time  of  maximum  winds  f  ~  35  hours,  the  Langmuir  cells  remain  coher¬ 
ent  and  are  readily  identified  in  the  images. 

The  impact  of  Langmuir  cells  in  the  high  wind  OBL  is  also  apparent  in  turbulence  statistics. 
Figure  4  compares  time  histories  of  vertical  velocity  statistics  from  three  different  LES  with 
varying  vortex  force,  viz.,  Lat  =  (oo,  0.5, 0.3).  Note  Lat  =  oo  is  a  simulation  with  no  vortex 
force  while  Lat  =  0.3  is  a  simulation  with  the  strongest  wave-current  interaction.  In  order  to 
construct  these  statistics,  we  first  compute  the  vertical  profiles  ( w2)(z,t )  and  ( w3)(z,t )  at  each 
time  step  of  the  LES;  here  (  )  denotes  a  horizontal  spatial  average  over  an  x  —  y  plane.  We  then 
sweep  through  these  profiles  to  locate  the  maximum  in  the  second-order  moment  ( w2)max  and 
minimum  in  the  third-order  moment  (w3)min.  The  temporal  variation  of  these  two  statistics  are 
given  in  figure  4.  Langmuir  cells,  generated  by  vortex  forces,  impact  the  vertical  velocity  and 
the  overall  turbulence  level  and  mixing  in  the  OBL  for  time  greater  than  10  hours.  Even  with 
relatively  weak  wave-current  interactions,  viz.,  Lat  =  0.5,  the  second  and  third  order  moments 
differ  from  the  simulation  with  no  vortex  force.  Grant  &  Belcher  (2009)  propose  Wl  =  u*/La^3 


as  the  velocity  scale  characteristic  of  the  Langmuir  turbulence  regime.  For  the  present  computa¬ 
tions  we  hold  the  turbulent  Langmuir  number  fixed  but  friction  velocity  u*  varies  considerably 
and  thus  the  simulations  sweep  through  a  wide  range  of  wl ■  Our  values  of  (w2)max  and  ( w3)min , 
obtained  with  time  varying  forcing,  are  broadly  consistent  to  those  obtained  by  Grant  &  Belcher 
(2009)  during  the  time  period  15  hours  <  t  <  40  hours.  One  of  the  future  tasks  is  to  examine  if 
the  velocity  scale  wl  collapses  vertical  profiles  of  important  OBL  turbulence  statistics  for  these 
high  wind  simulations. 

The  third-order  moment  (w3)  is  one  component  of  the  vertical  transport  term  T  =  — d{wq2)  /  dz 
which  appears  in  a  classical  turbulence  kinetic  energy  budget  ( q 2  is  the  total  turbulence  kinetic 
energy).  In  figure  5  we  show  vertical  profiles  of  Tw  =  —d(ww2)/dz  normalized  by  friction  ve¬ 
locity  and  boundary  layer  depth  h.  In  the  LES  with  no  wave  effects  this  transport  term  is  small 
with  little  vertical  structure  as  expected  for  a  shear  driven  boundary-layer  flow.  The  inclusion 
of  wave-current  interaction  in  the  simulations  generates  coherent  Langmuir  cells  that  alter  the 
turbulent  transport  in  striking  ways  (see  Grant  &  Belcher,  2009).  Significant  transport  of  ver¬ 
tical  variance  from  the  surface  layer  to  regions  deeper  in  the  OBL  is  induced  by  the  Langmuir 
cells.  This  enhances  the  overall  mixing  in  the  OBL.  For  example,  LES  solutions  show  that  the 
OBL  deepens  and  the  SST  is  cooled  in  the  presence  of  wave  effects  (not  shown).  However,  we 
find  that  most  of  the  change  in  these  two  statistics,  caused  by  wave  effects,  occurs  in  the  time 
period  before  the  time  of  maximum  winds. 


IMPACT  /APPLICATIONS 

The  LES  results  obtained  here  for  hurricane  driven  OBLs  can  be  used  to  guide  the  interpreta¬ 
tion  of  observations  collected  during  the  ITOP  program.  In  addition,  the  results  can  be  used  to 
test  simpler  1-D  parameterizations  of  the  ocean  mixed  layer  that  are  used  in  large  scale  models. 
Our  particular  interest  is  in  evaluating  and  improving  the  so-called  K-profile  parameterization 
(KPP)  that  is  routinely  used  in  the  Regional  Ocean  Modeling  System  (ROMS). 


TRANSITIONS  &  RELATED  PROJECTS 

The  present  work  has  links  to  the  ONR  DRI  on  High  Resolution  Air-Sea  Interaction  (HIRES) 
that  focuses  on  the  interaction  of  waves  and  turbulence  in  the  atmospheric  surface  layer.  The 
LES  model  being  developed  for  HIRES  is  also  being  used  in  the  present  work. 
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Figure  1:  Sketch  of  an  idealized  hurricane  vortex  used  for  setting  the  surface  conditions  in  LES.  The 
vortex  is  propagating  upward  with  speed  Vs  and  has  a  characteristic  length  scale  (radius  of  maximum 
winds)  A.  The  families  of  solid  and  dotted  circles  indicate  the  position  of  the  vortex  at  initial  time  t 
and  at  a  later  time  t  +  At,  respectively.  The  fixed  LES  domains  to  the  right  (R)  and  left  (L)  of  the 
vortex  center  (c)  feel  the  time  history  of  wind  speed  and  direction  along  the  dotted  vertical  lines. 
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Figure  2:  Temporal  variation  of  the  wind  and  thermal  forcing  applied  at  the  top  of  the  ocean 
boundary  layer.  The  LES  domain  in  figure  1  is  located  55  km  to  the  east  and  approximately  700 
km  north  of  the  hurricane  vortex  center.  The  vortex  propagates  vertically  at  5.5  m  s _1. 


Scalar  Flux  (w  m'2) 


Figure  3:  Contours  of  vertical  velocity  normalized  by  friction  velocity  u *  near  the  water  sur¬ 
face  at  different  times  in  the  LES  simulation  of  a  hurricane  driven  ocean  boundary  layer.  The 
elongated  downwelling  lines  (blue  contours)  mark  the  boundaries  between  neighboring  Langmuir 
cells.  Forward  looking  Y-junctions  indicate  mergers  between  neighboring  cells.  The  black  vec¬ 
tors  show  the  direction  of  the  surface  wind  applied  at  the  top  of  the  water  surface.  Notice  how 
the  direction  of  Langmuir  cells  mirrors  the  changing  rotation  of  the  wind.  In  these  simulations 
the  turbulent  Langmuir  number  Ldt  =  0.3. 


Figure  4:  Variation  of  the  maximum  vertical  velocity  variance  (upper  panel)  and  minimum 
third-order  moment  (lower  panel)  for  simulations  at  different  turbulent  Langmuir  numbers 
Lat=( oo,0.5,  0.3)  (blue,  green,  red)  curves  respectively.  The  time  of  maximum  winds  occurs  at 
about  t  =  36  hours  into  the  simulation.  The  variance  and  third-order  moment  are  normalized  by 
ul  and  ul,  respectively. 
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Figure  5:  Vertical  profile  of  the  (resolved)  vertical  transport  term  T=—d(ww2)/dz  which  ap¬ 
pears  in  the  turbulent  kinetic  energy  (TI\E)  equation.  The  vertical  transport  is  normalized 
using  the  ocean  boundary  layer  depth  h  and  friction  velocity  u * .  The  results  are  obtained  at 
t=35  hours  when  the  surface  wind  speed  S  ~  43.4  m  s-1.  The  winds  are  aligned  at  107  degrees 
from  the  horizontal.  The  turbulent  Langmuir  number  Lat=( oo,  0.5,  0.3)  (blue,  green,  red)  curves 
respectively.  Notice  in  the  case  with  strong  wave  effects  (red  curve)  turbulent  kinetic  energy  is 
transported  from  the  surface  layer  to  deeper  depths  in  the  boundary  layer  which  enhances  the 
entrainment  of  cool  water. 


